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Abstract: We report the finding of a scaling relation among the cosmic- web anisotropy pa- 
rameter A, the linear density rms fluctuation er(r) and the linear growth factor D(z). Using 
the tidal field derived from the Millennium Simulation on 512 3 grids at z = 0, 2, 5 and 127, 
we calculate the largest eigenvalues A of the local tidal tensor at each grid resolution and 
measure its distance-averaged two-point correlation function, £\, as a function of the cosines 
of polar angles cos 9 in the local principal axis frame. We show that £a is quite anisotropic, 
increasing toward the directions of minimal matter compression, and that the anisotropy 
of £a increases as the redshift z decreases and as the upper distance cutoff r c decreases. 
Fitting the numerical results to an analytic fitting model £a(cos#) oc (1 + Acos n 9)~ 1 , it 
is found that the best fit value of A, dubbed the cosmic-web anisotropy parameter, varies 
systematically with cr(r c ) and D(z), allowing us to determine the simple empiral scaling 
relation A(r c ,z) = 0.8 D°- 76 (z) a{r c ). 
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1. Introduction 

As confirmed by recent N-body simulations (l], |2|, the large-scale spatial distribution of cold 
dark matter exhibits an anisotropic web-like pattern, which is often dubbed the cosmic 
web. The cosmic web of the dark matter distribution found in N-body simulations is also 
consistent with the observed large-scale filamentary distribution of galaxies in the real uni- 
verse H ||] . According to the standard model of cosmic structure formation || , the cosmic 
web originates in primordial density perturbations that are sharpened by gravitational tidal 
fields. The web becomes more anisotropic as the tidal fields become stronger due to the 
nonlinear processes during the gravitational evolution. 

Various statistical tools have so far been proposed to describe the geometric proper- 
ties of the cosmic web. For example, the conventional Appoint statistics has been used to 
calculate the anisotropic spatial distribution of dark halos || . The Minkowski functionals 
have been employed to determine the morphological properties of the large-scale structures 
embedded in the cosmic web The A-dimensional skeleton approach has been found 
to be efficient in tracing the evolution of the cosmic web |, g. The Multiscale Morphol- 
ogy Filter method has been applied to the large-scale galaxy distribution to identify the 
anisotropic structures of the cosmic web [10]. A method utilizing the concept of the Lo- 
cal Dimension has been suggested for the local quantification of the shapes of the galaxy 



neighborhood in the cosmic web [11]. The ellipticity-ellipticity correlations of dark halos 
have been used to quantify the large-scale anisotropy of the cosmic web p^ |. Finally, tes- 
sellation techniques have been suggested to trace the evolution of the geometric structures 



in the universe |L3|, 14 1 



Very recently, it has been pointed out by Lee, Hahn and Porciani |15[|(LHP09, here- 
after) that, since the cosmic web is produced by the anisotropic compression of matter 



along the principal axes of the large-scale tidal fields, its anisotropic nature may be best 
quantified in the system of the principal axes of the tidal fields. This led them to suggest 
the anisotropic two-point correlations of the nonlinear traceless tidal field expressed in the 
principal axis frame, as a new statistical tool for the description of the cosmic web 

phenomenon. 

Analyzing numerical data from high-resolution N-body simulations, they have deter- 
mined £,\(r) an d found that it has a much larger correlation length (~ 20/i~ 1 Mpc) and a 
higher degree of anisotropy than the density field itself. Integrating £,\(r) over distance r 
from to a certain cut-off scale r c , and expressing it as a function of the angle between the 
major principal axes and the separation vectors, they noted that the nonlinear traceless 
tidal field has a much higher degree of anisotropy than the nonlinear density field. Inter- 
estingly, the results of LHC09 imply that the correlations of the nonlinear traceless tidal 
field may have a link to the initial conditions. 

Yet, the analysis of LHP09 was restricted to the present epoch and their results were 
obtained by setting the distance cutoff scale r c to the correlation length of the nonlinear 
tidal field, r c = 20/i~ 1 Mpc. In this paper, our goal is to explore how the two-point 
correlations of the nonlinear traceless tidal field expressed in the principal axis system 
vary with redshift and distance scale by analyzing numerical data from high-resolution 
cosmological simulations. Throughout this paper, we assume a flat ACDM cosmology. 

2. Physical Analysis 

2.1 Construction of the Tidal Fields 

For the construction of the tidal shear fields Tjj(x), we use the density contrast fields 6(x) 
constructed on 512 3 pixels from the Millennium Simulation by means of the count-in-cell 
method, at four different redshifts z = 0, 2, 5 and 127. The Millennium run [Q] followed 
the evolution of the trajectories of 10 10 dark matter particles, each of which has mass 
8.6 x 10 8 h^M Q , in a periodic box of linear size 500 h^Mpc, using a ACDM cosmology 
with the cosmological parameters f2 m = 0.25, = 0.75, h = 0.73, erg = 0.9 and n s = 1. 

The Fourier-transform of the density field, <5(k), is obtained through the Fast-Fourier- 
Transformation (FFT) method |jlq] . Then, the Fourier transform of the tidal shear field, 
Tjj(k), is calculated as 7ij(k) = kikj5()t) /k 2 . The inverse Fourier-transformation of Tjj(k) 
yields the tidal shear field in real space, T^-(x). The traceless tidal field is defined as 
Tjj(x) = Tjj(x) — 5{x)/3. At each pixel point we diagonalize iWx) to find its three eigen- 
values and the corresponding eigenvectors. The local principal axis frame has the three 
eigenvectors as basis vectors with polar axis in the direction of the eigenvector correspond- 
ing to the smallest eigenvalue. The eigenvector corresponding to the largest (smallest) 
eigenvalue is parallel to the direction of maximal (minimal) compression of local matter 
distribution. 

2.2 Correlations of the Traceless Tidal Fields 

The two-point correlation function of the largest eigenvalue A of the local tidal tensor at 



redshift z is defined as 



m 



6(r;2) = (A(x;z)-A(r + x;z)}. (2.1) 

Let us express the separation vector r in terms of the spherical polar coordinates in the 
system of the principal axes of the local tidal tensor, r = (r, cos 6, </>), with the polar 
axis aligned with the direction of minimum matter compression, i.e., the direction of the 
eigenvector corresponding to the smallest eigenvalue of Tjj(x). If the A-fleld was isotropic, 
then the two-point correlation (,\(r) would depend only on the separation distance r, but 
not on the polar and azimuthal angles 9 and <fi. The degree of the anisotropy of the cosmic 
web can be quantified by measuring how strongly £\ changes with 9 and <fi. As the 
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Figure 1: Anisotropic two point correlations of the lowest eigenvalues of the traceless tidal field 
in the tidal principal axis frame averaged over the distance up to r c as a function of the cosine of 
the angles between the directions to the neighbor points and the eigenvectors corresponding to the 
lowest eigenvalues at four different redshifts (z = 0, 2, 5 and 127 as dots, square dots, triangles 
and asterisks, respectively), for the four different values of the distance cut-off scale r c — 10, 15, 20 
and 25/i~ 1 Mpc. In each panel, the thin solid lines represent the fitting models. 



dependence of £,\(r) ° n the azimuthal angle <f> has been found to be rather weak 15], here 
we focus mainly on the polar-angle dependence of %)■ 

The anisotropic two-point correlation function of A as a function of the cosines of the 
polar angles can be obtained by averaging ° ver r an d 4> as 

Ucos9)=[ C dr[ d<££ A (r,cosM), (2.2) 
Jo Jo 

where r c denotes the upper distance cutoff. When we integrate £\(r) over r, we consider the 
distance cutoff scale r c greater than 10/i _1 Mpc since £a(cos#) shows fluctuating unstable 
behavior at r c less than 10/i _1 Mpc. It is expected that the degree of the anisotropy of 
£a(cos 9) depends on the distance cutoff r c and the redshift z, i.e., £a(cos 9) = £a(cos 9; r c , z). 

Using the data from the Millennium simulation, we numerically measure £a(cos 9; r c , z) 
for the cases of four different distance cutoff scales, r c = 10, 15, 20 and 25/i -1 Mpc, and at 
four different redshifts z = 0, 2, 5 and 127. For each pair of pixel points separated by r, 
we first compute the product of the largest eigenvalues of the local traceless tidal tensors 
T and determine the polar and azimuthal angles of r in the local principal axes of the tidal 
tensor. Then, we take the spatial average of it to calculate the correlation of the largest 
eigenvalue ^\(r;z). Finally, we determine the expression of the correlation as a function 
of the cosine of the polar angle by averaging £,\(r; z ) over the azimuthal angle 4> and over 
distance r from to r c . For the detailed explanation on how to measure £a(cos#) from 
numerical data, we refer the readers to LHP09 |15[ ]. 

Figure [l] plots the rescaled anisotropic two-point correlation of the traceless tidal fields, 
6 (cos (9) = £ A (cos0)/£\(O), for r c = 10, 15, 20 and 25/i" 1 Mpc (in the top-left, top-right, 
bottom-left and bottom-right panels, respectively) at z = 0, 2, 5 and 127 (dots, squares, 
triangles and open stars, respectively). In each panel, the solid lines correspond to the 
fitting models described in section 3.1. As it can be seen, £a(cos 9; r c , z) increases with 
cos 9 at all redshifts for all cases of r c , indicating that the A-field is anisotropic and more 
strongly correlated along the directions of minimum compression of dark matter in the 
local frame, which is consistent with the trend found by LHP09. 

The standard deviations on the numerical measurement of £a(cos#) are also calculated 
as statistical noises, which turn out to be as small as < 10 -5 for all cases. Although it would 
be desirable to estimate jackknife error bars that include not only the statistical noise but 
also the cosmic variance, it is difficult to estimate them here since dividing the simulation 
box would destroy the periodicity of the box needed for our Fourier-based calculation of 
the tidal field. Hence, we omit the errorbars in Fig. [l] since the calculated statistical errors 
are anyway invisibly small. 

3. Results 

3.1 Fitting Formula 

Noting that £a (cos 9;r c ) increases with cos 9, we employ the following fitting model for it 
©: 

£\(caB0.z:r„) = ; -. (3.1) 



where A and n are two adjustable parameters. Fitting Equation ( |3,l|) to the numerical 
results with the help of the x 2 - mm i mrza tion method, we determine the best-fit values 
of A and n for each case of r c at each redshift. For the calculation of x 2 , we set the 
values of all errors associated with the numerical measurements to unity rather than to the 
corresponding standard deviations of £\ [17], since the standard deviations are found to be 
extremely small in our case and do not account for cosmic variance, as discussed in section 
2.2. Table [l] lists the best-fit values obtained for A and n at z = for the four different 
values of r c , while in Table ^ we provide the best-fit values of A and n at four different 
redshifts, setting r c to 20/i~ 1 Mpc. 
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0.43 
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0.18 


2.21 
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0.34 


2.77 
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0.11 


2.26 
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0.28 


2.63 
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0.01 
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Table 1: Best-fit values of A and n at z 
for the four distance cut-off scales ry. 



Table 2: Best-fit values of A and n with 
r c = 20/i -1 Mpc at four different redshifts. 



Each panel of Fig. [T] shows the fitting models of £\ with the best-fit values of A and 
n as thin solid lines. It can be seen that Eq. (3.1) indeed provides excellent fits to the 
numerical results for all values of r c at all redshifts. In all cases the minimum values of x 2 
are formally found to be smaller than 10 -3 . The parameter A is dubbed the cosmic-web 
anisotropy parameter, as it measures the strength of the dependence of £\ on cos#. If the 
A-field was isotropic, A would be zero. The more anisotropic the cosmic web is, the higher 
the value of this parameter A. Regarding the other fitting parameter n, it also shows a 
monotonic decrease with r c , but no systematic change with redshift z. Therefore, we focus 
mainly on the cosmic web anisotropy parameter A from here on. 

It is worth comparing Eq. ( |3.1| ) with the fitting model given in LHP09, who defined 
the polar angle as the angle between r and the eigenvector corresponding to the largest 
eigenvalue. In this work, we define the polar angle as the angle between r and the eigen- 
vector corresponding to the smallest eigenvalue. The reason for using a different definition 
of the polar axis here is that £\ is found to exhibit more obvious anisotropic behavior when 
the polar axis is chosen to be parallel to the direction of the eigenvector corresponding to 
the smallest eigenvalue. In addition, we also vary the power of cos# in the fitting model, 
while in LHP09 this was kept fixed. 
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3.2 A Scaling Relation 

Noting a systematic trend in the variation of the cosmic- web anisotropy parameter A with r c 
and z, we investigate how A depends on the linear density rms fluctuation cr(r c ) smoothed 
on the distance cut-off scale r c , and on the linear growth factor D{z). Using various 
distance cut-off values equal to r c = 10, 12, 14, 16, 20 and 25 h Mpc at each redshift, we 
measure numerically £\(cos#;r c ) by using the data from the Millennium Simulation and 
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by determining the best-fit value of A through x 2 fitting of Eq. (|3.1| ) to the numerical data 
points in each case. 

We also analytically compute the linear density rms fluctuations cr(r c ) smoothed on 
the same distance cut-off scale r c as 

/(SO 
A 2 (k)W 2 (kr c )dlnk, (3.2) 
-00 

where W(kr c ) is the top-hat filter of scale radius r c and A 2 (/c) is the dimensionless linear 
matter power spectrum. The following analytic approximation for A 2 (fc) given by [|l^] is 
used: 

2 

[1 + 3.89<? + (16.1<?) 2 + (5.46) 3 + (6.71g) 4 ]~ 1/2 , (3.3) 

and n s is the spectral index of the primordial power 
spectrum. For the case that the key cosmological parameters (O m , h and as) are set at the 



A 2 (k) oc k ns+3 



m(l + 2.34q) 



2.34q 



where q = A;/[^ m /i 2 Mpc _1 ] [19 
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values used by the Millennium run, this formula turns out to agree sufficiently well with 
the result from the CMBFAST code [§(J. 

Figure || plots the best-fit values of A for the six different values of r c versus the 
linear rms density fluctuations <r(r c ) smoothed on the same scale, at the four different 
redshifts we considered. This plot reveals that the cosmic-web anisotropy parameter A 
is directly proportional to <r(r c ). The cosmic web anisotropy parameter A is shown to 
decrease monotonically with z, quantifying the evolution of the anisotropy of the traceless 
tidal fields. Note that this result is consistent with the qualitative explanation of the cosmic 
web theory ||. 

It is interesting to see in Fig. |2] that A[cr(r c )] at each redshift is well described by 
a straight line, and that the slope of the line decreases with z. Supposing A oc cr(r c ), 
we investigate how the proportionality factor between A and cr(r c ) (i.e., the slope of the 
straight line) changes with the linear growth factor D(z). For the evaluation of the linear 
growth factor D(z), we use the formula given in |2l|| . 

5 r°° 1 + 2/ 

D(2)a -fi,J [ Wl + 2 ) 3 + ^/ i ^ [SM1 + , 3 + iiAl 3 /2 - (3-4) 

We numerically calculate the ratio, A/a(r c ), at each redshift. Fitting A/a(r c ) to a 
power-law formula aD^(z) and determining the best-fit value of a and (3 at each redshift 
with the help of the x 2 statistics, we find that the ratio A/a(r c ) varies as 0.8 .D - 76 (z). 
Fig. H plots A(r c , z)/a(r c ) versus D(z). The numerical results of A(r Cl z)/a(r c ) at four 
redshifts are shown as circles while the fitting model 0.8 D 0,76 (z) is shown as a solid line. 
This result demonstrates that the numerically obtained ratio between A and <r(r c ) is indeed 
well fitted by 0.8D°- 76 (z). 

Finally, we can combine these findings into the following scaling relation among cosmic 
web anisotropy parameter A(r c ,z), linear density rms fluctuation o"(r c ) on the scale of r c , 
and linear growth factor D{z): 

A(r c ,z) = 0.8 D°- 76 (z)a(r c ). (3.5) 

This scaling relation quantifies how the cosmic web anisotropy parameter increases as 
the Universe evolves, revealing a link between the initial conditions and the anisotropic 
clustering of the nonlinear traceless tidal field expressed in the principal axis frame. 

4. Discussion 

To quantify the web-like pattern in the large-scale matter distribution, we have introduced 
the new concept of the cosmic web anisotropy parameter, which measures the degree of 
the anisotropy of the two-point correlations of the largest eigenvalues of the traceless tidal 
fields in the principal-axis frame. By analyzing the numerical data from the Millennium 
Simulation at different redshifts, we have empirically found a simple scaling relation among 
the cosmic web anisotropy parameter, the linear density rms fluctuations and the linear 
growth factor. Our results have allowed us to quantify how the anisotropy of the traceless 




tidal fields increases as the growth factor and the linear density rms fluctuations increase. 
This new scaling relation thus provides a measure for the growth of anisotropy in the 
evolving cosmic web of the Universe. 

The simple functional forms we found for the correlation and the scaling relation (eqs. 



1 3.1] and |3.5(| ) suggest that there might exist a simple explanation for these functional 
forms in linear perturbation theory. Yet, we have not been able to develop such a theory 
thus far, as the analytic treatment from first principles turns out to be extremely hard 
since we deal here with the nonlinear two-point correlations measured in the principal axis 
frame. 

A crucial implication of our result is that the anisotropy of the cosmic web may be pro- 
duced by three competing effects. The trace part of the tidal field and the cosmic expansion 
tend to make the matter distribution more isotropic, whereas the traceless part of the tidal 
field stretches the matter distribution and induces large-scale anisotropy. The competition 
among these three effects imprints the web-like pattern in the large-scale matter distribu- 
tion, just as the competition between gravity and radiation pressure has left imprints in the 



form of acoustic oscillations in the temperature map of the cosmic microwave background 
radiation [22]. 

Since the linear growth factor D(z) and the linear density rms fluctuations cr(r c ) are 
functions of the primary cosmological parameters such as the dark energy equation of 
state parameter w and the density parameter Cl m , our results also hint that the cosmic 
web anisotropy parameter may be a new probe of cosmology. Given the scaling relation 
(eq. |3.5| ), we expect that if the cosmic web anisotropy parameter A is measured for two 
different distance cut-off scales r c at the same redshift, it may be possible to put a constraint 
on the density parameter fi m by taking the ratio between the two values, removing the 
parameter degeneracy with the amplitude of the linear power spectrum a$. Similarly, by 
measuring the cosmic web anisotropy parameter at two different redshifts and taking the 
ratio between the two values, a constraint on the dark energy equation of state parameter 
w results. 

The success of using the cosmic web anisotropy parameter as a new cosmological probe, 
however, is contingent upon a couple of future tests. First of all, it needs to be examined 
whether or not the same scaling relation also holds for different cosmologies since we have 
derived it assuming a ACDM cosmology with a specific set of cosmological parameters. 
Secondly, what is readily measurable in practice is not the tidal field of the dark matter 
distribution but that of the galaxy distribution. It will hence be necessary to investigate 
how the bias between light and matter affects the scaling relation. We investigate this 
question in forthcoming work, and hope to report the results elsewhere in the near future. 
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